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We model s- and d-wave ceramic superconductors with a three-dimensional lattice of randomly 
distributed and vr Josephson junctions with finite self-inductance. The field and temperature 
dependences of the microwave absoption are obtained by solving the corresponding Langevin dy- 
namical equations. We find that at magnetic field H = the microwave absoption of the s-wave 
samples, when plotted against the field, has a minimum at any temperature. In the case of d-wave 
superconductors one has a peak at _H" = in the temperature region where the paramagnetic Meiss- 
ner effect is observable. These results agree with experiments. The dependence of the microwave 
absorption on the screening strength was found to be nontrivial due to the crossover from the weak 
to the strong screening regime. 
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PACS numbers: 75.40.Gb, 74.72.-h 



I. INTRODUCTION 

A fundamental property of superconductivity is the 
Meissner effect, i.e. the occurence of flux expulsion be- 
low the superconducting transition temperature and the 
resulting diamagnetic response to the external magnetic 
field. Contrary to this behavior a paramagnetic signal 
was observed in certain ceramic superconductors ut 
cooling in low enough fields (smaller than 0.1 mT) 
This effect is now referred to as the paramagnetic Meiss- 
ner effect (PME) or the Wohlleben effect. The nature 
of the unusual paramagnetic behavior may be related to 
the appearance q£ the spontaneous supercurrents (or of 
orbital moments)El. The latter occur due to the existence 
of TT-junctixans characterized by the negative |J|Qsephson 
couplingsH'Q. Furthermore, Sigrist and RiceHQ argued 
that the PME in the high- Tc superconductors is a con- 
sequence of the intrinsic unconveij^ional pairing symme- 
try of the HTCS of d^2^y2 typcE In fact, th&.PME is 
succesfully reproduced in a single loop modeH as well 
as in a model of interacting junction-loopttljj. The lat- 
ter model incorporates a network of Josephson junctions 
with a random concentration of 7r-junctions. The mag- 
netic screening is taken into account in both of the single 
and multi 7r-junction systems. 

The mechanism of the PME based on the d-wave sym- 
metry of the order parameter remains ambiguous be- 
cause it is not clear why this effect could not be ob- 
served in many ceramic materials. Furthermore, the 
paramagnetiG-resppase has been seen even in the conven- 
tional Nbll3ll3, Alllj^perconductors and the Nb-AlO^;- 
Nb tunnel junctionalj. In order to explain the PME 
in terms of conventional superconductivity one can em- 
ploy the idea of the flux compression inside of a sample. 
Such phenomenoD becomes possible in the presence of the 
inhoffljOgeneitieaij or of the sample boundarytj. Auletta 
et atll have also observed the PME in the model of spe- 
cial geometry involving only 0-junctions. In our opinion, 
the PME in this model is of the dynamical natMie but 
not the equillibrium effect as in the d-wave modelcl. Thus 



the intpinsic mechanism leading to the PME is still under 
debatem. 

One of the most valuable tools to distinguish between 
the s— and d-wave symmetry is the study of the m 
crowave absorption (MWA)ll3. In fact, Braunish et a, 
have found a nontrivial field dependence of the MWA in 
samples which display the PME. The MWA has a peak 
at 7? = 0, when plotted against H , whereas for s-wave 
superconductors it has a conventional minimum. Based 
on a hysteresis in the M — H space (M is a magnetiza- 
tion) Sigrist and Rice have shown that for the one-loop 
model the peak aX H = would be observed if the dimen- 
sionless self-inductance, L, exceeds some borderline value 
L* — 1. The relation between L and the inductance, L 
is as follows 



L = 



2nlc 



-L 



(1) 



where <I>o and Ic are the flux quantum and the critical 
current, respectively. j— . 

It should be noted that Dominguez et atB have qualita- 
tively reproduced the experimental findings for the MWA 
using the multi-loop model. Their results are, however, 
restricted to the two-dimensional system. More impor- 
tant, the question about the borderline value L* above 
which the nontrivial field dependence of the MWA may 
occur in the d-wave interacting loops model was not stud- 
ied. Also the role of temperature and of the screening 
have not been explored (the screening, for example, plays 
a key role in explaing experiments on the aging effect in 
ceramic superconductora23) . 

Since the underlying mechanism for the PME is still 
lacking, a careful study of MWA may shed some light 
on this problem. In this paper we study in detail the 
MWA in the three-dimensional system which is more rele- 
vant to experimental situations than the two-dimensional 
one. Integrating the corresponding Langevin equations 
we have made three new observations. First, contrary to 
the one-loop model the MWA in the system with ran- 
domly distributed 7r-junctions has a non-trivial field de- 



pendence for any value of L. In other words, L* — Q 
for the multi-loop model. Second, the peak at _ff = 
is found to disappear for T > T*, where T* is a boder- 
line temperature below which the PME is observable. T* 
was found to grow as the screening is lowered. The third 
observation is that for both s- and d-wave ceramics the 
MWA decreases with screening not monotonically but it 
has a minimum at L « 1. Such a behavior is related to 
a change in time and length scales when one goes from 
the weak to the strong screening limit. 



II. MODEL 

We neglect the charging effects, of the grains and con- 
sider the following Hamiltonianffl 

H = - ^ Jij cos(0, - 0j - Ai-j) + 

V 

^ i A, , A, = '-^ r Aif)dr , 



$. 



(2) 



where 6i is the phase of the condensate of the grain at 
the i-th site of a simple cubic lattice, A is the fluctuat- 
ing gauge potential at each link of the lattice, (pQ denotes 
the flux quantum, Jij denotes the Josephson coupling be- 
tween the i-th and j-th grains, L is the self-inductance 
of a loop (an elementary plaquette), while the mutual in- 
ductance between different loops is neglected. The first 
sum is taken over all nearest-neighbor pairs and the sec- 
ond sum is taken over all elementary plaquettes on the 
lattice. Fluctuating variables to be summed over are the 
phase variables, 9i, at each site and the gauge variables, 
Aij, at each link. $p is the total magnetic flux thread- 
ing through the p-th. plaquette, whereas ^p"^* is the flux 
due to an external magnetic field applied along the z- 
direction. 



(JjCXt 



HS if p is on the < xy > plane ,„, 

otherwise , ^ ' 



where S denotes the area of an elementary plaquette. 
For the d-wave superconductors we assume Jij to be an 
independent random variable taking the values J or —J 
with equal probability (± J or bimodal distribution), each 
representing and n junctions. In the case of s-wave 
ceramics Jij is always positive but distributed uniformly 
between and 2 J. 

It should be noted that model (2) is adequate to 
describe many dynamical phenomena plated to the 
PME such |-as the compesation effected, the aging 
phenomenoncJ, the effeC|L-pf applied electric fields, in the 
apparent critical currentcj and the ac resistivitycj. 



In order to study the MWA we have to calculate the 
linear response to an external electromagnetic field. Us- 
ing the relation between the MWA and the conductivity 
and the Kubo formulacJ one can show that this response 
is proportional to a voltage - voltage correlation function. 
Integrating over all of frequencies of the electromagnetic 
field we obtain the following expression for the MWA 



n = 



47r 
cnRT 



E 



<Vf > 



(4) 



where < V^ > is a mean value of the square of the voltage 
induced by the thermal noise on each junction, n is a light 
refraction coefficient and R is the normal resistance of the 
links. 

To calculate Vi we model the current flowing between 
two griSpiS with the resistively shunted junction (RSJ) 
model,ui23 which gives the following dynamical equations: 



h de^jn) 
2eR dt 



2e 



J^(n)sin6'^(n) 



2e£ 



A;[A+0^(n)-A+0,(n)]-7y,,(n,i) 



(5) 



Here we have redefined notation: the site of each grain is 
at position n = (n^,, rzy, n^) (i.e. i = n); the lattice direc- 
tions are /i = x, y, z; the link variables are between sites 
n and n -I- /i (i.e. link ij = link n, /i); and the plaquettes 
p are defined by the site n and the normal direction /i (i.e 
plaquette p = plaquette n, /i, for example the plaquette 
n, z is centered at position n -I- (x -I- y)/2). The forward 
difference operator A+6'y(n) = 0^(n + /i) — 9i,{n) and the 
backward operator A~0j,(n) = Ov{'^) ~ ^iy(n — /i). The 
Langevin noise current 77^ (n, t) has Gaussian correlations 



(r,^(n,i)v(n',t') 



2fcBT 
R 



S^,^'5n^n,S{t-t'). (6) 



The local voltage Vi is then given by 



K: 



d0^ 
dt 



(7) 



Eq. (5) describes the overdamped dynamics. We have 
tried to include the inertia (capacitive) terms but the re- 
sults do not change substantially and they are neglected. 

In what follows we will consider currents normalized 
by Ij — 2eJ/h, time by r = (Po/2ttIjR, voltages by RIj, 
inductance by 4)q/2ttIj and temperature by J /kB- Then 
the dimensionless MWA, Q, is defined as follows 



n 



cnR 

An 



n 



(8) 



III. RESULTS 

The system of differential equations (5) is inte- 
grated numerically by a second order Runge-Kutta- 
Helfand-Greenside algorithm for stochastic differential 



equationaHa. The time step depends on L and is equal to 
Ai = O.lrj and At = O.lrj x L for Z > 1 and X < 1, 
respectively. We conider the system size I — 8 (we have 
made some test runs for 1 — 12 and found that the finite 
size effects are not substantial). The temporal averages 
are taken over a time of 10^ tj after a transient time of 
25000tj. The free boundary conditions are implemented 
because the magnetizatiO|H_always vanishes for the peri- 
odic boundary conditionsBfl 
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FIG. 1. The field dependence of SI for s- (left panel) and 
d-wave (right panel) ceramic superconductors. We choose 
T — 0.2 and L = 0.1, 1 and 10. The results are averaged 
over 20 samples. 
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FIG. 2. The field dependence of the MWA for the d-wave 
samples. We took / = 8, L = 1 and T = 0.1, 0.5 and 1. The 
results are averaged over 10 -20 samples. 



Fig. 1 shows the field dependence of the MWA for 
r=0.2 and for various values of L. In the case of s-wave 
superconductors we have the standard minimum at H = 
for any value of inductance. It is also true for any T. As 
expected, fl ~ H^ at weak fields. For the d-wave samples 
fl has the unconveational peak at H — 0. Contrary to 
the one-loop modelQ such peak is seen not only for L > 1 
but also for L < I. It should be noted that the height 
of the peak is very small {{i^{H = 0) — i^min)/^min is 
of order of 10"'^). This is in qualitative agreement with 
experimental findingg3 that the peak should be low. 

In our model (2) the temperature dependence of the 
critical current is neglected. However, one can show that 
the dimensionless temperature T chosen in Fig. 1 and 
in all of the figures presented below corresponds to the 
relevant to experiments real temperature, T^. In fact, 
the critical current depends not only on temperature but 
also on conditions under which samples were prepared. 
The typical value of the cirtical current density for ce- 
ramic superconductors is ~ lO^A/m^ ( see, for exam- 
ple, Ref. |2^). Since the typical size of grains is about 
Ifim we have the critical current Ic ^ 10~^A. Using 
Tr = JT/kB = nicT/2ekB one obtains Tr/T - lOOif . 
Clearly, our dimensionless T cotsectly desribes the ex- 
perimental values of temperaturea. 
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FIG. 3. The inductance dependence of T* for d-wave su- 
perconductors. 

It is known that the random 7r-junction model {2^ dis- 
plays a phase transition to a so called chiral glass .l3 The 
frustration effect due to the random distribution of tt 
junctions leads to a glass state of quenched-in "chirali- 
ties" , which are local loop supercurrents circulating over 
grains and carrying a half-quantum of fiux. Evidence of 
this transition has been related to raeasurements of the 
nonlinear ac magnetic susceptibility.E3 The question we 
ask is if there is any correlation between the existence 
of the chiral glass phase and the anomalous behavior of 



il. As shown in Ref. |2^, the chiral glass disappears for 
L > Lc, where Lc — 5 — 7. On the other hand, the peak of 
il is obserable for any value of L. Therefore, there is no 
one-to-one correlation between the chiral glass and the 
nontrivial field dependence of the MWA of the ceramic 
superconductors . 

The field dependence of the MWA in the d-wave super- 
conductors for L = 1 and various values of T is shown 
in Fig. 2. At low T's the peak at H = shows up 
but it disappears at high temperatures. This is our main 
result. Such a observation was not reported in Ref. pj| . 
Qualitatively, above some borderline temperature, T* the 
frustration effect becomes less important and the d-wave 
system should behave like the s-wave one. 

Fig. 3 shows the dependence of T* on L. The ques- 
tion we ask now is why T* decreases with L. To answer 
this question we study the dependence of the roughness 
of the energy landscape on the screening. Since the num- 
ber of energy local minima should grossi-with the num- 
ber of grains (or of spins) exponentiallyE2l we restrict our 
calculations to small system sizes. We took I = 6 and 
search for the local minima by the annealing procedure 
at T > and then by the quenching at T = 0. The 
histogram, P{E), collected from local minima which are 
reached from 10000 starting configuations is shown in 
Fig. 4. Obviously, the energy landscape for L=l is more 
rugged compared to the L=10 case. 
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FIG. 4. The energy local minima histogram, P{E), for 
L = 1 and 10. We choose the system size I = 6. The energy 
bin used for collecting the histogram is equal to 0.002. The 
local minima are obtained from 10000 starting configurations 
by annealing at T > and quenching at T = 0. 

In order to characterize the roughness of the energy 
landscape we introduce the parameter S, 



V^Efr 



Im 



> 



< Elm. >2 



< Ei„i > 



(9) 



where Eim denotes the energy at local minima, < ... > 
means averaging over all minima studied. For the results 
presented in Fig. 4 we have 6 — 0.007 and 0.012 for 
L = 10 and 1, respectively. So the larger is screening the 
smaller roughness of the energy landscape. The differ- 
ence between s- and d-wave ceramics becomes, therefore, 
less and less pronounced as the screening increased and 
T* should go down with L. 

In order to understand the nature of T* we calculate 
the field cooled (FC) and zero field cooled (ZFC) mag- 
netization. In our model the magnetization is defined as 
follows 



M 



< 



1 



1(1 -ly 



p 



$0 



> 



HS 
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(10) 



where $p is the flux in a;y-plane, < ... > denotes the 
thermal and disorder average. In the FC runs, the tem- 
perature is lowered stepwise under a constant field. At 
each temperature, typically 10^ time steps are used for 
thermalization and 4 x 10^ steps for averaging. In the 
ZFC runs, the system is first quenched to a low temper- 
ature (T=0.05) in zero field and is thermalized during 
4 X 10^ steps. Then a static field is switched on and the 
temperature is increased stepwise under the same condi- 
tion as in the FC regime. 
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FIG. 5. (a)The temperature dependence of the FC and 
ZFC magnetization for the d-wave superconductors. L — 5. 
(b) The same as in (a) but for Cv (right-hand scale) and p2 
(left-hand scale). The latter was computed for the frequency 
of the ac electric field u; — 0.001 and its magnitude Jo = 0.1 
(see Ref. 23 for details). The results are averaged over 20 - 
40 samples. 

Fig. 5a shows the FC and ZFC magnetization for L = 5 
at a finite magnetic field / ~ HS/(I)q = 0.1. We can see 
that T* is the temperature below which one has an onset 
of positive magnetization, i.e. the paramagnetic Meiss- 



ner effect starts to be observed. The irreversibility point 
occurs at temperatures lower than T* , and its position 
is dependent on the heating or cooling rate. We identify 
T* to correspond to the onset of chiral short-range or- 
der where pair-loops begin to appear locally. The main 
conclusion here is that the anomalous field dependence 
of the MWA is strongly correlated with the occurence of 
the PME. On the other hand, as was shown in Ref. |[ the 
PME may appear for any value of screening, we conclude 
that for the multi-loop model the borderline value L* for 
the anomalous behavior of MWA is equal to 0. 

Fig. 5b shows the temperature dependence of the spe- 
cific heat, C^, and the nonlinear ac resistivity p2 for 
L = 5. Ci, is defined as porportional to the energy fluc- 
tuations, Cv = {{5EY)/T'^. The definition of p2 is given 
in Ref. ESl Clearly, T* coincides with the peak of C„ and 
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FIG. 6. The inductance dependence of the MWA at i/ = 
for T = 0.2 (squares), 0.5 (hexagons) and 0.8 (circles). The 
system size I = 8. The results are averaged over 10 samples. 

We now study the depedence of Vt on the screening. 
Our results are shown in Fig. 6 for H = but the quali- 
tative behavior is also valid for H ^{). There is no appre- 
ciable difference between s- and d-wave cases. For a fixed 
value of L, Q.{H = 0) depends on T very weakly. The de- 
pendence on screening is more pronounced. From naive 
arguments the MWA should decreases with L because the 
screening would prevent the absorption in the bulk. Fig. 
6 shows, however, that fi has a minimum at L = 1. The 
anomalous dependence of MWA on L may be understood 
in the following way. The static asd-dynamic properties 
of the Joscphson arrays are shownEl'E^ to change qualita- 
tively if the inductance varies from L < 1 to L > 1. The 
attractive vortex-vortex interaction in the weak screen- 
ing regime becomes repulsive in the opposite limit. The 
qualitative change in the dynamic response is related to 
the change of length and time scales. Since the magnetic 



screening length goes as A ^ L~^/^ (see Ref. 25), for 
L < 1 A is larger than the grain size (lattice spacing in 
the cubic network) while for L > \ \ becomes smaller 
than the grain size. For L < \ the relaxation time for 
the field is smaller than the relaxation time for the phases 
whereas the opposite happens for L > 1. The decrease 
of the phase relaxation time compared to the field one 
should, therefore, increase the MWA for L > 1. 



IV. CONCLUSION 

In conclusion, experimental results of Braunish et al.u 
for the MWA can be reproduced by the XY-like model 
for the d-wave superconductor. Although the peak of 
il is found to be small its study is useful for elucidat- 
ing the symmetry of the superconducting order parame- 
ter. Within the multi-loop model the anomalous behav- 
ior should be observable for any value of inductance if 
T < T* . At high temperatures there is no qualitative 
difference between the s- and d-wave systems. The de- 
pendence of the MWA on the screening strength is found 
to be not monotonic due to the crossover from the weak 
to the strong screening regime. It would be very inter- 
esting to verify this prediction experimentally. 
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